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We consider a gas of trapped Cooper-paired fermionic atoms which are manipulated by laser 
light. The laser induces a transition from an internal state with large negative scattering length 
(superfluid) to one with weaker interactions (normal gas). We show that the process can be used 
to detect the presence of the superconducting order parameter. Also, we propose a direct way of 
measuring the size of the gap in the trap. The efficiency and feasibility of this probing method is 
investigated in detail in different physical situations. 
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I. INTRODUCTION 

Recent experiments on cooling and trapping Fermionic 
atoms have opened np new opportunities for study- 
ing fundamental quantum statistical and many-body 
physics. Trapped Fermionic '*''K atoms were cooled down 
to temperatures where the Fermi degeneracy sets in [Q. 
Two lithium isotopes were trapped simultaneously in 
an magneto-optical trap in 1^ and optical trapping of 
fermionic lithium has been achieved as well . The rich- 
ness of the internal energy structure of the atoms and the 
possibility to accurately and efficiently manipulate these 
energy states by laser light and magnetic fields allows ex- 
cellent control of these gases. Furthermore, atomic gases 
are dilute and weakly interacting thus offering the ideal 
tool for developing and experimentally testing theories of 
many-body quantum physics. 

The degenerate Fermi gas is expected to show many 
interesting phenomena in its thermodynamics , excita- 
tion spectrum l5-|^, coUisional dynamics |^ and scatter- 
ing of ligh t p] , p^ ,pd|. A major goal is to observe the pre- 
dicted 13 1 BCS-transition for Fermionic atoms - this 
would compare to the experimental realization of atomic 
Bose-Einstein condensates It is still, however, an 

open question how to observe the BCS-transition, be- 
cause the value of the superconducting order parameter 
(gap) is expected to be small and the existence of the or- 
der parameter does not significantly change the density 
profile and other bulk properties of the gas. 

There are several proposals for measuring the super- 
conducting order parameter. Off-resonant light scatter- 
ing as a probe was proposed in [|l5|,^. Superfluidity 
is predicted to effect both the spectral and spatial dis- 
tribution of the scattered light. In |17|, superfluidity 
was found to increase the optical lineshift and lincwidth. 
Also non-optical phenomena, such as collective and sin- 
gle particle excitations, have been proposed to be used 
for observing the BCS-transition Jl8|-^. Probing by a 
magnetic field was considered in [ESL 

The use of on-resonant light as a probe for the order pa- 



rameter was proposed in . The basic idea is to transfer 
atoms from one internal (hyperfine) state for which the 
atoms are Cooper paired to another state for which the 
interatomic interaction is not strong enough to lead to a 
BCS state. This effectively creates a superconducting - 
normal state interface across which the atomic popula- 
tion can move. There is a conceptual analogy to electron 
tunneling from a superconducting metal to a normal one 
which is used to measure the gap and the density of states 
for Cooper-paired electrons . The tuneability of the 
interaction strengths which is required for this scheme 
is obtained by the use of magnetic fields. This allows 
to manipulate the scattering lengths between atoms in 
different internal states, see e.g. the recent experimental 
results concerning optically trapped fermionic Lithium 
atoms 1^ and the theory predictions for '^'^K |26|] . 

The basic idea in the proposal is that the absorp- 
tion peak is shifted and becomes asymmetric because of 
the existence of the gap - the laser has to provide en- 
ergy for breaking the Cooper pairs in order to transfer 
atoms from the paired state to the unpaired one. This 
behaviour is, however, strongly influenced by the spe- 
cific physical situation. In this paper we investigate in 
detail how the choice of the chemical potentials for the 
superfiuid and the normal state, and the choice of the 
interaction strengths and laser profiles affect the absorp- 
tion. We also compare the results in the cases of a ho- 
mogeneous system and a trapped gas. In section we 
introduce the considered system. The linear response of 
the gas for a light probe is derived in section |ll^, b oth 
for a homogeneous and a trapped gas. In section IV the 



various parameters that affect the observed absorption 
are discussed. In section ^ exact numerical results are 
presented for the limit when the laser beam profile can 
be considered a constant. A beam profile with non-zero 
intensity only in a small volume in the middle of the trap 
is considered in section This case is very interesting 
because then only the center of the trap is probed. Thus 
the order parameter seen by the laser is almost constant 
and indeed, we find a remarkable agreement with the re- 
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suits predicted for a homogeneous system. We finally 



summarize the results in section VII 



II. 



THE SUPERCONDUCTOR - NORMAL STATE 
INTERFACE 




FIG. 1. Probing of the gap in a gas of attractively interact- 
ing cold Fermionic atoms, a) Laser excitation with the cou- 
pling fl and the detuning A transfers a Cooper paired atom 
from the internal state \g) to the state |e). b) The other atom 
in the initial Cooper pair becomes an excitation in the BCS 
state, therefore the laser has to provide also the additional 
gap energy Ac- In this picture the Fermi levels Hg and fie for 
the internal states have been chosen to be different from each 
other but they could also equal. 

We consider atoms with three internal states available, 
say |e), \g), and \g'). They are chosen so that the inter- 
action between atoms in states \g) and \g') is relatively 
strong and their chemical potentials are nearly equal so 
that \g) and \g') can be assumed to be Cooper-paired. 
All other interactions are small enough and/or the chem- 
ical potentials of the corresponding states are different 
enough in order to assume the |e) atoms to be in a nor- 
mal state [Q. The laser frequency is chosen to transfer 
population between |e) and \g), but is not in resonance 
with any transition which could move population away 
from the state \g'). 

For small intensities, the laser-interaction can be 
treated as a perturbation, the unperturbed states being 
the normal and the superconducting state. The transfer 
of atoms from \g) to |e) is then analogous to tunneling 
of electrons from a normal metal to a superconductor in- 
duced by an external voltage, which can be used as a 
method to probe the gap and the density of states of 
the superconductor . In our case the tunneling is be- 
tween two internal states rather than two spatial regions, 
this resembles the idea of internal Josephson-oscillations 
in two-component Bose-Einstein condensates . Fig. 1 
illustrates the basic idea. The observable carrying essen- 
tial information about the superconducting state is the 
change in the population of the state |e), we call this the 
current /. 

We define a two-component Fermion field 



(1) 



where ipe and ipg fulfill standard Fermionic commutation 
relations. The fields ipe/g can be expanded using some 



basis functions (e.g. plane waves or trap wave functions) 
and corresponding creation and annihilation operators: 



V'e/g(x) = c^'^^(/'^''^(x). The annihilation and creation 

t 



operators fulfil {cf ,c^^} v -^o- , j- — uij. 
two components of the field, corresponding to the inter- 
nal states |e) and \g), are coupled by a laser. This can 
be either a direct excitation or a Raman process; we de- 
note the atomic energy level difference by Ua {h = 1), the 
laser frequency lol and the wave vector - in the case 
of a Raman process these are effective quantities. In the 
rotating wave approximation the Hamiltonian reads 

Here A — Wa — ujl is the (effective) detuning, and r2(x) 
contains the spatial dependence of the laser field mul- 
tiplied with the (effective) Rabi frequency. The parts 
He and Hggi contain terms which depend only on ipe or 
Ipg, Ipg', respectively. Possible spatial inhomogeneity, e.g. 
from the trap potential, is also included in and Hggi . 



Oand {cl^^lcfj^S,.^ 



The 



III. THE CURRENT 

The observable carrying essential information about 
the superconducting state is the rate of change in popu- 
lation of the |e} state. We may call it, after the electron 
tunneling analogy, the current 

I{t) = -{Ne), 

where = J d'^xipl{x.)ipe{^)- The current I{t) is calcu- 
lated considering the tunneling part of the Hamiltonian, 

HT^H-{He+ Hgg, + A/2(7Ve - Ng)) 

as a perturbation; the current / becomes the first order 
response to the external perturbation caused by the laser. 
We calculate it both in the homogeneous case and in the 
case of harmonic confinement. The calculations are done 
in the grand canonical ensemble, therefore the chemical 
potentials fj,g and /ig are introduced. Also the detuning 
A acts like a difference in chemical potentials, thus it be- 
comes useful to define an effective quantity of the form 
A = /ie — /ig + A = A/i + A. In the derivation we assume 
finite temperature, but most of the results will only be 
quoted for T=0. 



A. Homogeneous case 

The assumption of spatial homogeneity is appropriate 
when the atoms are confined in a trap potential which 
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changes very little compared to characteristic quantities 
of the system, such as the coherence length and the size 
of the Cooper pairs. In the present context, this assump- 
tion is also valid when the laser profile is chosen so that 
it only probes the middle of the trap where the order 
parameter is nearly constant in space - this will be dis- 



cussed in more detail in section Vl 



In the homogeneous case the Fermion fields V'e/g can 
be expanded into plane waves. The Hamiltonian becomes 



H = H, + H„ 



k 



kl 



(3) 



r2(x) = rie^^^'^. The term \Tki\'^ now produces a delta- 
function enforcing momentum conservation. Note that 
this is very different from the assumption of a constant 
transfer matrix (X^fei l^fc'l^ — ^ I^I^Sfcz) made in the 
standard calculation for tunneling of electrons over a su- 
perconductor - normal metal surface The final re- 
sult becomes (assuming for simplicity that the tempera- 
ture T = 0) 



-7rl7V(A)6'(-A 
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where k is given by the following energy conservation 
condition: 



where 
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We calculate the current 



(4) 



treating Ht as a perturbation: terms of higher order 
than are neglected. Because we are interested in the 
current between the superconducting and normal states, 
correlations of the form {clclcgCg) (and h.c.) are omitted 
since they correspond to tunneling of pairs (Josephson 
current). The current can be written 



dt9{t){e-'^'{[AHQ),Am-e^^\[A{Q),AHt)])) 



where 



and 



m = EkiTkici\t)cf{t), 

where K = H - Ht - fieNe - HgNg. The 



two terms in the above equation have the form of retarded 
and advanced Green's functions. These are evaluated us- 
ing Matsubara Green's functions techniques, which leads 
to 
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-nF{e+A)]Ag{k,e+A)Ae{l,e), 



where np are the Fermi distribution functions and Ag/^ 
are the spectral functions for the superconducting and 
normal states. We use the standard expressions [p5|] 



A,{l,e) = 2TrS{e-Ci) 



and 



Ag{k, e + A) = 2Tr[ulS{e + A - ujk) + vlS{e + A + ujk)]- 

Here = Ei — fi,, and u^., and ujk sue given by the 
Bogoliubov transformation. Here we consider for sim- 
plicity the term proportional to u^; is analogous. 
The laser field is chosen to be a running wave, that is 



-A- 



k — kr 



tk-k. = 0- 



= VCl + Aj, and 



P(A) 



V 

2^ 



A 



is the density of states which appears when the summa- 
tion over momenta is changed into an integration over 
energies. 

The laser momentum kj^ can be very small compared 
to the momentum of the atoms, especially in the case of 
a Raman process. By setting kL — Q the result becomes, 
including now terms proportional to both vf, and u\: 



Ah 



2AmA)^, 



(5) 



where ± are for A > and A < 0, respectively. The 
term with the — sign corresponds to A < i.e. current 
from \g) to |e), and the positive term to current from 
|e) to \g). To understand the results in terms of physics, 
let us first consider the case of equal chemical potentials 
A^ = 0: 



A2 



^nQ'^piAM-A- 



eiA-Ac)]. (6) 



In order to transfer one atom from the state \g) to |e) 
the laser has to break a Cooper pair. The minimum en- 
ergy required for this is the gap energy A^, therefore the 
current does not flow before the laser detuning provides 
this energy - this is expressed by the first step function 
in (||). As |A| increases further, the current will decrease 
quadratically. This is because the case |A| = Aq corre- 
sponds to the transfer of particles with p = pp, whereas 
larger |A| means larger momenta, and there are simply 
fewer Cooper-pairs away from the Fermi surface. This 
behaviour is very different from the electron tunneling 
where the current grows as ^/{eVy~~Kh (the volt- 
age eV corresponds to the detuning A in our case) be- 
cause all momentum states are coupled to each other. 
The second step function in (H) corresponds to tunneling 
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into the superconductor. In this case one has to provide 
extra energy because a single particle tunneling into a 
superconductor becomes a quasi-particle excitation with 
the minimum energy given by the gap energy. 

When the chemical potentials are not equal the situa- 
tion is more complicated, but the basic features are the 
same: i) threshold for the onset of the current determined 
by the gap energy and difference in chemical potentials, 
and ii) further decay of the current because the density 
of the states that can fulfill energy and momentum con- 
servation decreases. 



B. Harmonic confinement 

In the case of harmonic confinement the spatial depen- 
dence of the current is non-trivial. We define the total 
current as 



where A^jg are defined 



lit) = - / d3^(7Ve(x)), 



(7) 



where 



iVe(x) = A^e(x)] = l[r!*(x)^t(x)^^(x) - ff.c.]. (8) 

No expansion in the plane waves or other basis functions 
is made at this point and the first order perturbation 
calculation leads to a result with explicitly spatially de- 
pendent correlation functions: 

/ = 2/m[X,et(-A)], (9) 

where 

X^et(-A) = i f dte-'^*9{t) f (fx 

J —OO J 

X / d3a;'([At(x,0),^(x',i)]), (10) 



and 



A(x,t) = f]*(x)V4(x)V'e(x). 



(11) 



The retarded Green's function Xret(— A) is calculated us- 
ing Matsubara techniques. For performing the Matsub- 
ara summations, the spatially dependent Green's func- 
tions for the normal and the superfluid state are ex- 
panded using the trap wave functions (/)„(x), and the 
Bogoliubov-de Gennes (BdG) eigenfunctions u„(x) and 
z;„(x), respectively. This leads to the result 



2^ 



(fx / d^x'n*{x)n{x) 



[ie(x,x',e)G^,,(x',x,6 + A) 
+G^,t(x,x',e- A)ig(x',x, e) 



i,/,(x,x',e) = z(G?/f(x,x',e) - Gl2i^y,e)) 



and 



M„(xO<(x) ^ <(xOi;„(x) 



fx x' - V M^Wni^') 



(13) 
(14) 



As mentioned, u„(x), Un(x) and ujn are given by the 
BdG equations [p8[ . In taking the imaginary part of 
the expression (|12[) we first collect together all spatially 
dependent terms, which gives real factors of the form 
I J d^xQ{x)un{x)(j)m{^)\'^ (wc also use the fact that the 
trap wave functions (/)m(x) are real). Imaginary parts of 
the remaining terms give spectral functions in the usual 
way. The derivation leads to 



I = -2ttY^ y d^2:l](x)u„(x)(/)„(x) 



(i^a;r2(xX(x)(/)™(x) 



[npi-i^n) - f^F(?m)]<5(Cm + A -|- w„). 



(15) 



(12) 



Note that this form docs not lead to a simple step- 
function type behaviour like in the homogeneous case. 
Due to the non-orthogonality of the trap and the BdG 
wavefunctions transitions between many quantum num- 
bers are allowed and the total current is the sum of all 
these. 

In the following sections we use the result (|l^) to in- 
vestigate the feasibility of the method as a probe in dif- 
ferent physical situations. The eigenfunctions and val- 
ues w„(x) and LUn are calculated numerically from the 
BdG equations using the pseudo-potential method pre- 
sented in 1^^. We assume for simplicity that the trap 
has spherical symmetry. The quantum numbers n, m in 
Eq.(|l5|) then become 77, 1, m with I, ni being the usual an- 
gular momentum quantum numbers. The quasi particle 
(QP) energies will only depend on r/, I. The method of 
solving the BdG equations is described in detail in [p9| . 

Note that the above derivation assumes that the BdG 
equations are solved exactly. One can, however, also use 
the local density approximation ||30| where the chemi- 
cal potential of the superconducting state is assumed to 
have a parametric dependence on position. In this case 
one should repeat the above derivation by assuming that 
the superfluid state Hamiltonian Hggi is parametrically 
depended on x: Hgg,. This leads to a result identical to 
( [Tsl ) except that now u„ , w„ are actually plane waves but 
they as well as LUn have a parametric dependence on x 
via ^(x). We may denote the current for a chosen x as 
The total current is then the average of for all x. 
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IV. THE EFFECT OF CHEMICAL POTENTIAL, 
LASER PROFILE AND TRAPPING 

The behaviour of the current is strongly influenced by 
the choice of the physical parameters. This allows a con- 
venient way to optimize the probing scheme as well as 
to investigate interesting physics such as the influence of 
the harmonic confinement — this is our twofold aim in 
the following. 

Although for instance the laser profile can be chosen 
by will, there are a few basic restrictions in choosing 
the chemical potentials and the interaction strengths be- 
tween the atoms in different internal states. In order for 
a gas of Fermionic atoms in two internal states to form 
Cooper pairs, the interatomic interaction should be large 
enough and the chemical potentials corresponding to the 
two states should be very close to each other Q. This 
is the condition we assume for \g) and \g'). The state 
|e) is always assumed to have either negligible interac- 
tion with the states \g) and \g'), or considerably lower 
chemical potential (smaller number of particles) . It turns 
out that not only the pairing effects the light absorption 
but also normal interactions described by the Hartree- 
field are crucial. This is illuminated by a comparison of 
these two cases of |e) having negligible scattering length 
or small chemical potential. 

In real experiments, the gas is trapped in a harmonic 
potential. We have derived results also for the homoge- 
neous case, these can be used in the case of very large 
traps and they also give a simple intuitive picture of the 
basic physics in this system. We will show that indeed the 
trapping has considerable effect on the results. On the 
other hand, an effectively homogeneous situation can be 
achieved by choosing the laser to probe only the center of 
the trap. As will be shown, this avoids certain problems 
arising from the non-homogeneous potential and can give 
a very clear signature of the superconducting state. 



V. CONSTANT BEAM PROFILE 

In this section we present exact numerical results in 
the limit when the laser beam intensity can be well ap- 
proximated by a constant. That is, we take r2(x) = in 
Eq. (|l5|) . As the typical effective wavelength for the laser 
is much larger than the extend of the trap for the relevant 
energies, such an approximation should be good as long 
as the laser amplitude is constant over the whole extend 
of the cloud. We present results for two situations: the 
case when there are no atoms initially in the state |e) and 
the case when the chemical potentials for |e} and \g) are 
the same. 



A. No [e) atoms initially 

In the case considered in this subsection, we assume to 
have initially only a gas of interacting \g) and \g') atoms. 
The laser beam then induces transitions to the hyperfine 
level |e). To obtain clear results, we will assume that the 
e) atoms see the same Hartree field as the \g) atoms. 
This means that geg + geg' — ggg' since the |e) atoms 
see the Hartree field from both the \g) and \g') atoms. 
Here, gggi = ATragg'/m denotes the interaction strength 
between the two hyperfine states \g) and \g') and like- 
wise for geg and geg>- The parameter agg> is the usual 
s— wave scattering length for scattering between \g) and 
\g') atoms. Experimentally, this situation could possi- 
bly be achieved by manipulating an external magnetic 
field thereby tuning the effective low energy interaction 
between the relevant hyperfine states to an appropriate 
value. In Fig. ^, we show a typical example of the cur- 
rent /(A) — —{Ne). We have used parameters such that 
gggi — —I'lgUj, fig = 3 1 . 5 w and the temperature T = 
with /ig denoting the chemical potential for the \g) and 
\g') atoms. Here Iho = (mw)""'^/^ is the harmonic oscilla- 
tor length. For ^Li atoms with a^g' =-1140A this cor- 
responds to approx. 1.6 x 10^ atoms trapped in the states 
\g) and \g') with a trapping frequency of 820Hz yielding 
a critical temperature Tc ^ llOnK for the BCS transi- 
tion. We have added an imaginary part T — O.lhu! to 
the quasiparticle energies such that the 7r(5(x)-functions 
in Eq.(|l|) become Lorentzians T/2{x^ + r^/4). 

Normal-normal current The dashed curve in Fig. ^ 
depicts the current when the \g) and \g') atoms are in 
the normal phase. Since the |e) atoms see the same 
Hartree-field as the \g) atoms, the spatial part of the QP 
wavefunctions is the same for the two hyperfine states. 
Thus, the overlap integrals in Eq.(|l5|) simply become 
(5„m-functions and since the quasiparticle spectrum is 
the same for the two hyperfine states, we obtain from 
Eq.®: 

^(A) = -2^. Al^- (16) 

where Ng is the number of \g) atoms trapped. 

Superconductor-normal current The solid curve Fig. ^ 
depicts the current when \g) and \g') are in the superfluid 
phase. We see that the maximum of the current is dis- 
placed from A = and that the shape of the current pro- 
file is asymmetric. Both effects are quite straightforward 
to understand. The asymmetry reflects the fact that the 
current now is given by a sum of Lorentzians centered at 
different frequencies since the QP spectra for \g) and |e) 
are different and the overlap integral in Eq. (^|) does not 
give a simple selection rule. The shift in the center of 
the peak to negative A is due to the fact that in order 
to induce a transition from a low lying QP-state in the 
superfluid one needs to break a Cooper pair. This re- 
quires an additional energy given by the pairing energy 
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of the QP-state. As a fraction ~ T^Tp of the particles 
participates in the pairing and they on average have the 
pairing energy one can estimate the order of mag- 
nitude of the shift in the center of the peak away from 
its normal phase value A = to be 0{T^ /Tp). Here, 
Tc is the critical temperature for the BCS transition and 
ksTp = fj,g is the Fermi temperature for the \g) atoms. 
For the present parameters, we have fcsTc ~ 2.8^ giving 
T^/Tp ^ 0.25 in qualitative agreement with the results 
depicted in Fig. (||). We have performed numerical cal- 
culations for several values ggg> and ^g, and we find the 
general behavior as described in the present example. In 
all the tested cases, the current peak for the superfluid 
phase is shifted away to negative values of A and the 
shape of the peak is asymmetric as opposed to the sim- 
ple Lorentzian shape for the normal phase. The shift of 
the peak is of the order 0(T^/Tp). 

In order to enhance the effect of the pairing on /(A) 
even further, one could initially also trap some |e) atoms. 
As long as fig — fie » Ac, the |e) atoms will not 
Cooper pair with the \g) or \g') atoms even though g^g or 
geg> < By having the lower QP states for the |e) 
atoms filled, there will only be transitions between the 
QP states around the Fermi energy. Since these states 
are influenced the strongest by the pairing, the effect of 
the superfluidity on /(A) will be even stronger than for 
the parameters relevant for Fig. ^. This is illustrated by 
Fig. ^, where we plot /(A) for the same parameters as 
above apart from that now /ig = 21. bTiu. We see that 
both the asymmetry and the shift in the peak in the su- 
perfluid phase as compared to the normal phase, is more 
pronounced than in Fig. (^). This is simply because the 
transitions deep below the Fermi level which are essen- 
tially inert to the effects of superfluidity, are blocked by 
filling up the levels for the |e) atoms up to the energy 
E = 21.bhu. However, it might be somewhat more dif- 
ficult to achieve this situation experimentally, as it re- 
quires the initial trapping of 3 (instead of 2) hyperfine 
states with a rather good control of the populations in 
each state. 

One should note that it is important that the Hartree 
field seen by the \g) and |e) atoms is approximately the 
same. Otherwise, the wavefunctions and the spectra for 
the \g) and |e) atoms will be different even when the \g) 
and l^') atoms are in the normal phase. The overlap in- 
tegrals will then not give simple selection rules and there 
will be a contribution from many Lorentzians centered in 
general away from A = 0. Consequently, /(A) will not 
be given by the simple formula in Eq. (p^. This is illus- 
trated in Fig. ^, where we plot the current /(A) for the 
same parameters as given above (with no |e) atoms ini- 
tially), apart from we now have geg + geg' — 0.9 x ggg'. As 
expected, the current profile is shifted away from A = 
and is asymmetric, even when the \g) and \g') atoms are 
in the normal phase. The shift to negative frequencies 
is easy to understand: The attractive mean (Hartree) 



field seen by the |e) atoms is smaller than the attractive 
field seen by the \g) atoms since geg + geg' = 0.9 x gggi. 
Therefore, the trap states for the |e) atoms in general 
have a slightly higher energy than for the \g) atoms, and 
the normal phase current is shifted to negative A. Fig- 
ure ^ demonstrates that the pairing field still causes a 
general shift of /(A) to negative frequencies and intro- 
duces further asymmetry since the pairing energy still 
needs to be broken to generate a current from the su- 
perfluid phase. This effect is readily visible since the 
Hartree fields seen by the \g) and |e) atoms are approx- 
imately the same for the parameters chosen. However, 
if the difference in the Hartree fields becomes too large, 
the spread in the signal is determined by this difference 
and any additional effects coming from the pairing field 
are correspondingly obscured. In general, to be able to 
detect the presence of superfluidity using the scheme de- 
scribed in this section, one should have AW ^ Aq where 
AW = \ggg' — geg — geg'\p denotes the difference in the 
Hartree fields with p being the average density of the \g) 
atoms, and Aq is the average gap. 

We conclude that if geg + geg' = ggg' to a very good 
approximation so that the difference in the Hartree fields 
seen by the \g) and the |e) is negligible, the effect of 
superfluidity on the current /(A) should be straightfor- 
ward to observe. The current in the normal phase is a 
simple Lorentzian centered around A = 0, whereas in 
the superfluid phase it is asymmetric and shifted away 
from A = 0. Furthermore, the shift in the center of the 
peak provides an estimate of Te if Tp is known. Both 
the asymmetry and the shift away from A = should be 
easily observable indications of the presence of superflu- 
idity. The effect is further enhanced if one initially traps 
|e) atoms keeping pg — fXe 3> Aq. In general, the scheme 
described in this section requires that the difference in 
the Hartree fields seen by the \g) and the |e) is smaller 
than the average pairing field in order to obtain a visible 
effect of the superfluidity on the current. 

B. Equal chemical potentials 

We now consider the case of /Zg — /Xg, where there is 
initially many atoms trapped in the state |e). Hence, to 
avoid the otherwise interesting possibility of the |e) atoms 
participating in the pairing, we assume that geg — geg' = 
0. In Fig. we plot the current /(A) for ggg> — —I'l^uj, 
fie = fJ-g ~ 31.5a;, T=0, and geg = geg' = 0. We have 
taken F — O.Ilu. As can be seen, there are several peaks in 
/(A) both when \g) and \g') are in the normal phase and 
when they are in the superfluid phase. In both cases, the 
peaks simply correspond to individual QP energy bands 
overlapping. 

To see this, we plot in Fig. |^ the corresponding lowest 
QP energies for the |e), \g), and \g') atoms in both phases 
as a function of their angular momentum I. As we are 
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in the Bogoliubov picture, all QP energies are positive 
and measured relative to the chemical potential. In the 
normal phase, negative QP energies ^„ = e„ — ^ < are 
represented as positive QP energies En = |^n| with hole 
character. In Fig. (H), we label a hole-state by o whereas 
a particle state is indicated by a x. In the superfluid 
phase, the QP's are in general a superposition of a hole 
and a particle which we label by +. 

We see that when the Hartree field is attractive, a nor- 
mal phase energy band with a downward curvature in 
Fig. ^ is a hole band whereas a normal phase energy band 
with an upward curvature is a particle band. The reason 
is that particle states with lower angular momentum I 
has a lower energy for an attractive interaction (g < 0) 
than states with a higher I since the wavefunction overlap 
with the Hartree field decreases with increasing I [Q. 

The QP bands for the je) atoms are flat as they are 
the simple unperturbed harmonic oscillator states with 
energies En = \{n + 3/2)uj — He\- The lowest E = 
band corresponds to the harmonic oscillator states at 
the chemical potential {n — 30) with angular momen- 
tum 1 = 0,2..., 30. The interpretation of the spectra is 
described in detail in |^,^^. There is a one to one corre- 
spondence between the peaks in Fig. |5| and the QP bands 
depicted in Fig. |^. For instance, the broad peak centered 
around A = uj when the \g) and \g') atoms are in the 
normal phase corresponds to transitions from the half- 
filled QP band at = for the |e) atoms into the empty 
(particle) band with 0.5 < E/oj < 1.8 for Z = 0, 2 . . . 30 
for the \g) atoms in the normal phase. 

We note that the peaks for the superfluid phase are 
sharper than the peaks for the normal phase. This is be- 
cause the lowest QP bands for the \g) atoms are almost 
degenerate as a function of I in the superfluid phase as 
can be seen from Fig. ^. These states are strongly in- 
fluenced by the pairing field, which "pushes" them away 
from the center of the trap. They are concentrated in 
the region between where the pairing field and the trap- 
ping potential are significant ||l^,^, and they thus do 
not feel the Hartree field. Therefore, their dependence 
on the quantum number I is much weaker than in the 
normal phase. 

We have performed calculations for a number of differ- 
ent parameters and we have observed the general behav- 
ior of /(A) as described above. We conclude that when 
we have /ie — ^J'g and geg = geg> = 0, the presence of 
superfluidity is somewhat harder to detect as compared 
to the situation described in sec. VA. This is because 



the Hartree field tends to obscure any additional effect 
coming for the pairing. The current /(A) has in general 
many peaks corresponding to the energetic overlap be- 
tween individual QP bands for the |e) and the \g) atoms 
both for the normal and the superfluid phases. The effect 
of the superfluidity is to make the peaks sharper than in 
the normal phase since the pairing field tends restore the 
degeneracy in the QP energies with respect to the angu- 



lar momentum. One could therefore possibly detect the 
onset of superfluidity as a sharpening of the peaks. 



VI. PROBING THE CENTER OF THE TRAP 
EFFECTIVELY HOMOGENEOUS SYSTEM 



We now assume that the laser intensity 51 (x) is large 
in the center of the cloud and that it decreases quickly 
as a function of the distance from the center of the trap. 
This situation can be experimentally achieved by using 
a Raman transition scheme with two perpendicular laser 
beams crossing each other at the center of the cloud. The 
profile of each laser beam should be narrow on the length 
scale of the trapped cloud. Since the laser beams then 
effectively probe only atoms in the center of the cloud 
where the Hartree and pairing fields are approximately 
constant, we would expect the observed signal to be well 
described by the results for a homogeneous system as 
given by Eq.(|^). From section 



HI A, we conclude that, 



for a homogeneous system, the optimal way of detect- 
ing the presence of the pairing field is to have fig ~ /ie 
(see Eq.(^). In this limit, all low lying transitions far 
away from the Fermi energy and thus very little affected 
by the pairing, are Fermi blocked. The transitions con- 
tributing to /(A) are all close to the Fermi level and 
hence strongly infiuenced by the presence of the pairing 
field. We will therefore concentrate on the case where the 
effective chemical potential in the center of the trap is the 
same for \g) and |e). As we will see, this case opens up 
the interesting possibility of directly measuring the size 
of the gap in the center of the cloud. 



We use the same set of parameters as in sec. VB for 
the \g) and \g') atoms. But now we assume that the 
intensity profile for the beam can be well approximated 
by a sphere of constant intensity for r < ro and zero in- 
tensity for r > vq with r denoting the distance to the 
center of the trap. That is, we take 0(r) = r28(ro — r) 
in Eq.(15) with ro ~ 2lho- /_,From Fig. 0, we see that 
AG(r) ~ 6cj and W{r) ~ 15w for r < 2lho. Thus, the 
effective local potential for the \g) atoms is /ig ~ 46.5a;. 
We therefore take fig = 46.5a; and g^g = geg' = in order 
to have the same effective local chemical potential for \g) 
and Ie) in the center of the trap. The result is shown in 
Fig. @. Since the system is approximately homogeneous 
for r < 2lho, one could expect the current profile to be 
well described by Eq.(|^), with fi^ = l^-g = 46.5a; and 
Ag = 6a;. We therefore also plot the result predicted by 
Eq.(^. Here we take fie = fJ-g — 46.5a; and Aq = 6a; and 
we have normalized /(A) to a volume oi V = 47rro/3. 
We have taken a rather large value of F = la; such that 
the discrete nature of the trap spectrum is washed out. 
Note that a finite imaginary part F to the QP energies 
corresponds to convoluting Eq.(|^) with a Lorentzian of 
width r. 
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As can be seen, there is a good agreement between 
the exact numerical result and the prediction based on 
Eq.(^. Especially, the current is zero for — Gti) < A < 6lu 
as predicted by Eq.(^), since one either needs to break a 
Cooper pair with pairing energy ~ Aq (r — 0) to produce 
a current into |e) (=> A < — 6w) or one has to create a 
QP with energy minimum ~ Ac^r = 0) A ^ 6uj) to 
generate a current into \g). The peaks in the numerical 
result reflect, of course, the discreteness of the trap levels. 
These peaks are quite large, as there is only w 500 par- 
ticles trapped in the region r < 2lho for the parameters 
given above. Clearly, these peaks cannot be reproduced 
by the homogeneous treatment given by Eq.(|^) which 
however reproduces the general shape of the current pro- 
file well. If we had chosen a larger system the individual 
peaks would be more numerous and smaller on the scale 
of the gap Aoi^r = 0) and the agreement between the 
homogeneous approximation and the exact result would 
probably be even better. 

We conclude that by concentrating the beam inten- 
sity to the center of the cloud where the gas can be to a 
good approximation regarded as homogeneous, the cur- 
rent profile /(A) is well described by the results presented 
in section III A. An important result is that by adjusting 



the parameters such that the local chemical potentials 
are the same [fig — Wg{r = 0) ~ /ie — We{r = 0)], one 
should be able to measure directly the size of the gap in 
the center of the cloud. It is simply given by the thresh- 
old detuning energy below which the observed current 
should be zero: for |A| < Aq the current /(A) ~ 0. 



the normal state. Furthermore, the scattering length be- 
tween the normal state atoms and the Cooper paired ones 
is assumed to be about the half of the scattering length 
between the two Cooper paired ones - this causes all the 
atoms to see the same Hartree field. In this physical sit- 
uation the effect of the BCS-state is particularly simple 
and clear: the maximum in the current of population as a 
function of the detuning is shifted and the peak becomes 
asymmetric. Although this would probably be the opti- 
mal choice, other initial conditions and probe parameters 
lead to clear signatures of the BCS-state as well. 

To avoid the problems arising from the non- 
homogeneous trapping potential we propose to probe 
only the middle of the trap. The order parameter is effec- 
tively homogeneous in the middle and the wave- functions 
of the in-gap excitations are located away from the cen- 
ter. In practice this kind of probing can be done by 
using two orthogonal Raman beams which intersect only 
in the middle of the trap. We have shown that indeed 
this leads to a result very similar to the one in the homo- 
geneous: the maximum of the current is shifted exactly 
by the amount of the gap energy. This allows a direct 
measurement of the gap energy. 
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VII. CONCLUSIONS 



The observation of the predicted BCS-state in gases of 
trapped atomic Fermions poses a double challenge: the 
order parameter is small thus a sensitive probe has to 
be found, furthermore, the trapping potential leads to 
appearance of in-gap low-energy excitations which may 
make it difficult to resolve the gap energy. In this paper 
we have presented a method which is based on the trans- 
fer of atomic population over a superconductor - normal 
state interface. This interface is effectively created by us- 
ing a laser to couple internal states with large and small 
scattering lengths. The population transfer requires to 
break a Cooper pair and the extra energy for this is pro- 
vided by the laser detuning. The change in the atomic 
population as a function of the laser detuning thus gives 
information about the gap energy. We have derived the 
current of population both in the case of a trapped gas 
and a homogeneous system, and investigated the feasi- 
bility of the method in different physical regimes. 

We found that, in the case of a constant laser pro- 
file, the clearest signatures of the BCS-state are observed 
when it is assumed that there are initially no atoms in 
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FIG. 2. The current / — —{Ne) as a function of the 
laser detuning for = 31.5, = (in units of uj), and 

geg + Seg' ~ 9gg' ■ 
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FIG. 3. The current I = -(Ne)- The sohd/dashed lines 
are for the \g} and \g') atoms in the superfluid/normal phase 
for fig — 31.5, /le = 21.5 (trap units), and geg + deg' = dgg'- 
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